Parity- violating asymmetry of W bosons 
produced in p-p collisions 



Xun Chen a , Yajun Mao a , Bo-Qiang Ma bA * 

^School of Physics, Peking University, Beijing 100871, China 
h CCAST (World Laboratory), P.O. Box 8730, Beijing 100080, China 



Abstract 

The parity-violating asymmetry is an ideal tool to study the quark helicity dis- 
tribution in the proton. We study the parity-violating asymmetry of bosons 
produced by longitudinally polarized p-p collision in RHIC, based on predictions of 
quark distributions of the proton in the SU(6) quark-spectator-diquark model and 
a perturbative QCD based counting rule analysis. We find that the two models give 
nearly equal asymmetry for W + but that for W~ quite different. Therefore future 
experiments on such quantity can help to clarify different predictions of the value 
Ad(x)/d(x) at x — ► 1 in the proton. 
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1 Introduction 



The helicity distribution is one of the three fundamental quark distributions 
of the nucleon, and it has been explored for more than two decades in deep 
inelastic scattering progresses. Deep inelastic scattering (DIS) experiments 
of polarized electrons and muons from polarized nucleons have confirmed us 
that the total quark-plus-antiquark helicity sum is only about 1/4 to 1/3 of 
the proton spin [1,2,3,4,5]. However, inclusive DIS experiments cannot provide 
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information on the polarized quark and antiquark density separately. Semi- 
inclusive DIS measurements are one of the approaches to achieve a separation 
of quark and antiquark helicity densities [6,7]. But these experiments depend 
on the details of the fragmentation process, so their accuracy is limited. In 
these experiments, spin asymmetry is one of the most important quantity to 
study the quark helicity distributions. The Jefferson Lab Hall A Collaboration 
has released a precision result on the neutron spin asymmetry [8,9] at large 
x, which provides an examination on some known quark distribution models. 
The disadvantage of semi-inclusive DIS process can be avoided by using the 
production of W bosons at RHIC [10,11]. 

About ten years ago, Soffer and his collaborators pointed out that the po- 
larization of u, d, u, and d in the proton will be measured directly and pre- 
cisely using the maximum parity violation of W bosons in ud — > W + and 
du — > W~ at high Q 2 [14,15]. At leading order, the single spin asymmetry, or 
the parity- violating asymmetry, of W ± will approach to the ratio Aq(x)/q(x) 
when the absolute value of the rapidity of W, yw, is large. Some works have 
been done to calculate this asymmetry with different quark distribution mod- 
els [12,13,14,15,16,17,18,19,20]. The purpose of this work is to show that the 
parity-violating asymmetry of W~ bosons in longitudinally polarized p-p pro- 
cess at RHIC can provide a crucial test of different predictions on Ad(x) / d(x) 
from a pQCD based analysis and the SU(6) quark-spectator-diquark models, 
because it is very sensitive to the value of Ad(x) / d(x) at x — > 1 in the proton. 

The parity- violating asymmetry is defined as the difference of the cross section 
of the left-handed (— ) beams and that of the right-handed (+) beams, divided 
by the sum of them, or using the number density of the beams instead and 
normalized by the polarization: 

AL ~ a + + a_~ P X N>, + N>_- (1) 

The production of W is dominated by u, d, u, and d contributions, with some 
contamination of s, c, s, and c contributions, mostly through quark mixing. 
Therefore W production is an ideal tool to study the spin-flavor structure of 
the nucleon. The expected sensitivities at RHIC are also given in Ref. [10]. 
PHENIX and STAR both expect to get the enough W + (W~) data for analysis. 

At RHIC, the W ± will be produced in sub-processes ud — > W + and du — * W~ 
by colliding beams of protons spinning alternately left- and right-handed. 

For Standard Model, the leading order production of W + s, is partly illustrated 
in Fig.l [10]. Due to that the production of left-handed W bosons violate 
parity, the asymmetry can be constructed from polarized beams with two 
directions . As Fig. 1 shows, the u quark only comes from the polarized proton, 
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Proton helicity=" +" 



Proton helicity="-" 





Fig. 1. A possible production of W + in p-p collision, at lowest order. The u quark 
comes from the polarized proton and the d quark comes from the unpolarized proton. 

the asymmetry should be expressed as: 



A w+ = u_{x 1 )d{x 2 ) -u + (xi)d(x 2 ) = Aw(xi) 
L uZ(x 1 )d(x 2 ) + u+(xi)J(x 2 ) u(xi) 



(2) 



In the above equation, superscripts refer to the quark helicity and subscripts 
refer to that of proton. q(x n ) is short for q(x n , Q 2 ), where q = u,d and n — 1,2, 
Q 2 = Myy and Mw is the mass of W bosons. 

For the situation that the d quark comes from the polarized proton, the asym- 
metry should be simply written as: 



w + _ dt{xi)u(x 2 ) - d+(xi)u(x 2 ) _ Arf(ari) 
dt(xi)u(x 2 ) + dX{xi)u(x 2 ) d(xi) 



(3) 



In general, the asymmetry is a superposition of the two cases: 
Au(xi)d(x 2 ) — Ad(xi)u(x 2 ) 



aT = 



u(x 1 )d(x 2 ) + d(x 1 )u(x 2 ) 



(4) 



The asymmetry of W can be got easily by interchanging u and d: 
Ad(xi)u(x 2 ) — Au(xi)d(x 2 ) 



A w- 

L 



d(xi)u(x 2 ) + u(x 1 )d(x 2 ) 



(5) 



Higher-order corrections change the asymmetries only a little [16,17]. 

The kinematics of W production and Drell-Yan production of lepton pairs 
are the same. The momentum fraction carried by the quarks and antiquarks, 
x n {n = 1,2), can be determined by the rapidity of W, yw, 



Mi 



Xi 
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x 2 = 



W - 



yw 



(6) 



where y/s is the center-of-mass energy. This picture is valid only for the pre- 
dominant production of Ws at pr = 0. In reality, the W has pr, which comes 
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Fig. 2. Predictions on g^/Ff of quark-diquark model [23] and pQCD based analy- 
sis [26] for Q 2 = 4GeV 2 , and are compared to the experimental data from HER- 
MES [5] and E155 [3]. 

from the higher level contributions such as gu — > W + d and ud — > W + g, or 
the primordial pt for the initial parton. But we must point out that when 
y w > 1, AY + « Au(x)/u(x), and ~ Ad(x)/d(x). So if we take p T = 
in calculation, the results are also valuable at large yw 

The experimental difficulty is that the W is observed through its leptonic 
decay W — > lis, and only the charged lepton is observed. The connection of 
the rapidity of observed leptons and that of Ws has been described in detail 
in Ref. [10]. And the expected cross section of Ws has also been estimated in 
that article. 



2 Quark Distributions 

From above equations, we can calculate the asymmetries of W^s if the quark 
helicity distributions are given. Two known models, the quark-diquark model [21,22,23] 
and a pQCD based counting rule analysis [24,25,26,27,28], can provide pre- 
dictions of the quark distributions. Both models have their advantages and 
have played important roles in the investigation of various nucleon structure 
functions. However, there are still some unknowns concerning the sea con- 
tent of the nucleon and the large x behaviors of valence quark. For example, 
there are still some uncertainties concerning the flavor decomposition of the 
quark helicity distributions at large x, especially for the less dominant d va- 
lence quark of the proton. There are two different theoretical predictions of 
the ratio Ad(x)/d(x) at x — > 1: the pQCD base counting rule analysis [26] 
predicts Ad(x)/d(x) — > 1 whereas the SU(6) quark-diquark model [23] pre- 
dicts Ad(x)/d(x) — > —1/3. Experimental data of recent years show that both 
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Fig. 3. Predictions on g^fF™ of quark-diquark model [23] and pQCD based analy- 
sis [26] for Q 2 = 4 GeV 2 , and are compared to experimental data from JLab [8,9] 
and HERMES [4]. 
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Fig. 4. Predictions on the flavor decomposition (An + Au)/(u + u) and 
(Ad + Ad)/(d + d) of quark-diquark model [23] and pQCD based analysis [26] 
for Q 2 = 4 GeV 2 , and are compared to experimental data from JLab [8,9] and 
HERMES [7]. 

of these two models are good in the predictions of g\, and of the Au(x)/u(x) 
extracted from it. But results of Jefferson Lab Hall A Collaboration show that 
the pQCD prediction on A\ deviates from the data [8,9]. Figures 2,3 and 4 
show the predictions of these two models and the comparison with experi- 
mental data. It seems that the pQCD predictions fit the data good at small 
x < 0.3 whereas the SU(6) quark-diquark model works good at large x > 0.3. 

We can expect that the different predictions on Ad(x)/d(x) of the two models 
may lead to prominent differences in calculating the parity- violating asymme- 
try of W~ , due to the case that AY ~ Ad(x)/d(x) when yw is large. 
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Detailed constructions of the quark helicity distributions in the two models 
can be found in Ref. [29]. 



To make realistic predictions of measurable quantities, we need also to take 
into the sea quark contribution in the two model constructions. In the quark- 
diquark case, this can be achieved by adopting one set of unpolarized quark dis- 
tribution parametrization as input, and then use theoretical relations to con- 
nect the quark helicity distributions with the unpolarized distributions [23,30]: 

Au v (x)=[u v (x) - ^d v (x)]W s (x) - ^d v (x)W v (x), 

1 ^ 
Ad v (x)=--d v (x)W v {x). 



Ws(x) and Wy(x) are the Melosh-Wigner rotation factors [31] for spectator 
scalar and vector diquarks. We use the valence quark momentum distributions 
u v (x) and d v (x) from quark distribution parametrization, but with Ws(x) and 
WV(a;) from the model calculation [30]. In this way we can take into account 
the sea contribution for the unpolarized quark distribution from the input 
parametrization. The energy scale dependence and the Q 2 evolution behaviors 
of the quark distributions can be reflected by the explicitly Q 2 dependence 
of the input quark distribution parametrization. This can provide a more 
reliable prediction for the magnitude and shape of a physical quantity than 
directly from the model calculation. CTEQ DIS [32] (CTEQ5 set 2) quark 
distributions are used for parametrization to investigate the evolution of the 
helicity distribution. From Fig. 5 we found that the ratios of Au ( x )+ A ^) anc [ 

A rf(xj+J(zT' > changed little with the increase of Q 2 in this model. 

For the pQCD based analysis, we take the same consideration as the above and 
make the following connection to relate the pQCD model quark distributions 
with the parametrization: 



u^ CD (x)=u^(x), (8) 

rf pQCD ( ^ ) = ^M M para (x)) (9) 

Au^ CB (x) = ^^u^(x), (10) 

ArfP QCD (x) = A4jM M para (a;)) (n) 
ul n (x) 

where the superscripts "th" means the pure theoretical calculation in the 

pQCD analysis, and "para" means the input from parametrization. The su- 
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Fig. 5. The Q 2 evolution of the helicity distributions. Prediction of the decompo- 
sition (Ait + An) /{u + u) and (Ad + Ac?) /{d + d) of quark-diquark model (thick 
curves) and pQCD based analysis (thin curves) for Q 2 = 1 GeV 2 (solid curves) and 
50 GeV 2 (dashed curves). 



perscript "pQCD" means that the new quark distributions keep exactly the 
same flavor and spin structure as those in the pQCD analysis, but with detailed 
x-dependent behaviors more close to the realistic situation. This is equivalent 
to using a factor, -uF ara \x) / 'u^ 1 , to adjust each pure theoretically calculated 
quantity to a more realistic pQCD model quantity. In this way we can take 
into account the sea contribution by using the sea quark distributions from 
parametrization, while still keep the pQCD model behaviors of the valence 
quark distributions. We also use CTEQ DIS quark distribution functions as 
input to investigate the Q 2 evolution of helicity distributions. The results are 
drawn in Fig. 5. We found the ratios changed little with the increase of Q 2 in 
this model as those in SU(6) quark-diquark model. 

Thus we have two set of quark distributions of Aq(x) and q(x), which keeps 
the same valence behaviors as in the quark diquark model and pQCD based 
analysis prediction. In the SU(6) model, the sea quarks are unpolarized, and 
only valence quarks are polarized. Both of these two models do not give the 
polarized distribution of the sea. HERMES data [33] show that the absolute 
value of polarized sea quark distributions is much smaller than the absolute 
value of polarized valence quark distributions when x > 0.1. So it does not 
matter for large x and yw if we omit terms including Aq(x) in the calculation. 
The Q 2 evolution behavior of the decomposition indicates that the asymmetry 
changes little when calculated with distribution functions at different Q 2 , thus 
the results and conclusion of this work will not influenced qualitatively by the 
Q 2 evolution effect. 
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Fig. 6. The value of AY (yw) m different models. The solid curves correspond to 
the SU(6) quark-spectator-diquark model, and the dashed curves corresponds to 
the pQCD based analysis, for yfs = 200 GeV. 




Fig. 7. The value of AY (yw) m different models. The solid curves correspond to 
the SU(6) quark-spectator-diquark model, and the dashed curves correspond to the 
pQCD based analysis, for -y/i = 500 GeV. 

3 Numerical Calculations 

With the definition of and the quark distributions, we can perform nu- 
merical calculations using different models. For the unpolarized quark distri- 
butions, we use the CTEQ LO parametrization [32] (CTEQ5 set 3) as input 
for both the SU(6) quark-spectator-diquark model and pQCD based analysis 
distribution functions described above. Q = M w = 80.419 GeV is taken to 
get quark distributions from the CTEQ LO parametrization, where Mw is the 
mass of W bosons. 

The polarized p-p collisions at RHIC will take place at center-of-mass energies 
oi yfs = 200 ~ 500 GeV. So we calculate A% at y/s = 200 GeV and ^ = 
500 GeV. Although the production of W at 200 GeV is not feasible with the 
designed luminosity at RHIC, we provide the calculation as a reference to 
show the model dependent predictions at different collision energies. Results 
of AY (yw) are shown in Fig. 6 and Fig. 7. 

It is clear that AY + are nearly equal in these two models, but AY show great 
difference between the two models at both the two energy scales. 

As we mentioned before, AY ~ Ad(x)/d(x) when yw is large, which means 
that x\ is large. We may see the difference more clearly by taking x\ as the 
independent variable. To examine this, we calculate AY again. Results are 
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Fig. 8. The value of A^(xi) in different models. The solid curves correspond to 
the SU(6) quark-spectator-diquark model, and the dashed curves correspond to the 
pQCD based analysis, for yfs = 200 GeV. 




Fig. 9. The value of (xi) in different models. The solid curves correspond to 
the SU(6) quark-spectator-diquark model, and the dashed curves correspond to the 
pQCD based analysis, for = 500 GeV. 

shown in Fig. 8, where */s = 200 GeV, and Fig. 9, where */s = 500 GeV. 

The results convince us that the difference of in these two models comes 
from the different predictions on Ad(x)/d(x) of the proton. We expect that 
experiments at RHIC can give an examination. 



4 Summary 



The parity-violating asymmetry in the production of W bosons in longitudi- 
nally polarized p-p process can be measured at RHIC. It is sensitive to quark 
distributions of the proton. A pQCD based analysis and the SU(6) quark- 
spectator-diquark model give obvious different results on the parity-violating 
asymmetry of W~, which is sensitive to the value of Ad(x)/d(x), a basic quan- 
tity which is difficulty to be measured in other processes. Therefore RHIC 
experiments results can give a powerful examination on different predictions 
and enrich our knowledge on the spin structure of the proton. 
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